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It was recently claimed that on d-dimensional small- world 
networks with a density p of shortcuts, the typical separation 
s(p) ~ p~ 1//d between shortcut-ends is a characteristic length 
for shortest-paths (M. E. J. Newman and D. J. Watts, "Scal- 
ing and perco lation in the small-world network model' 



cond- 



mat/990441^)- This contradicts an earlier argument suggest- 



ing that no finite characteristic length can be defined for bilo- 
cal observables on these systems (M. Argollo de Menezes, 
C. Moukarzel and T. J . P. Penna, "First -order transition in 
small- world networks" , cont-mat/990342t ) . We show analyti- 



cally, and confirm by numerical simulation, that shortest-path 
lengths £(r) behave as £(r) ~ r for r < r c and as £(r) ~ r c 
for r > r c , where r is the Euclidean separation between two 
points and r c (p,L) ~ log(L d p). This shows that the 

mean separation s between shortcut-ends is not a relevant 
length-scale for shortest-paths. The true characteristic length 
r c (p, L) diverges with system size L no matter the value of p. 
Therefore no finite characteristic length can be defined for 
small- world networks in the thermodynamic limit. 

PACS numbers: 05.10.-a, 05.40.-a, 05.50.+q, 87.18.Sn 



The recent observation jj| that a non-extensive 
amount of long-range bonds, or "shortcuts", on a reg- 
ular ci-dimensional lattice suffices to drastically shorten 
shortest-path distances, has triggered appreciable inter- 
est in these so-called "small- world" networks || ||. 

Small-world net- 
works consist of a regular d-dimensional lattice of linear 
size L to which L d p long-ranged bonds have been added, 
connecting randomly chosen pairs of sites. Several recent 
studies have concentrated on the behavior of the average 
shortest-path distance I — l/-/V 2 ^<r,-> < % >) where 
£ij is the minimum number of links that must be tra- 
versed to join sites i and j, and <> means average over 
disorder realizations. On a regular lattice one has £ ~ L, 
while on a random graph of L d sites, £ ~ log(i). On 
small-world networks, which interpolate between these 
two limits, one finds Jl|,|,| 1 ~ L if L d p « 1, but 
I ~ logL if L d p >> 1. Thus on a large system a small 
density of shortcuts is enough to have shortest-path dis- 
tances which are characteristic of random graphs. 

It has been suggested || that for any fixed density p of 
long-ranged bonds, a crossover-size L*(p) exists, above 
which the average shortest-path distance I increases only 
logarithmically with L. Renormalization Group argu- 
ments H as well as numerical measurements in- 
dicate that this crossover length diverges as p — > as 
L*(p) ~ p' 1/d . While some authors see this result as ev- 
idence of a continuous phase transition at p = feUfl] , it 



has also been argued M that no finite correlation length 
can be defined on these systems, and thus that the small- 
world transition is first-order. 

Clearly, the normalized shortest-path distance £(p) = 
£/L undergoes a discontinuity at p = in the thermody- 
namic limit. Thus, the "small-world" transition is tech- 
nically a first-order, or discontinuous, transition 0. A 
subtler and still controversial point is whether it 

is possible to identify a correlation length diverging 
at p — 0, i.e. if the small-world transition is associated 
with some sort of critical behavior. This would mean that 
p = is a first-order critical pomt(FOCP) JTlQ . As dis- 
cussed by Fisher and Berker a FOCP is a critical 
point where the thermal eigen-exponent attains its max- 
imum possible value y t = d, thus allowing the coexistence 
of two different critical phases. Consequently at a FOCP 
one has v — l/yt = 1/d, and some finite-size corrections 
will be dictated by this exponent. 

A FOCP appears for example in one-dimensional per- 
colation JTo| , where the order parameter (the den- 
sity of the spanning cluster) is discontinuous at p c — 1 , 
while the correlation length £ diverges as p — » 1 as 
£ ~ (1 — p)^ 1 ■ Finite-size effects as e.g. the width of 
the critical region are of order L _1 . Because of this, 
is non-zero even for p 1, ifL<£~(l — p)^ 1 . Another 
example of a FOCP is the Ising model in one dimension, 
where the magnetization density < m > (T) is discon- 
tinuous at T — while the correlation length £ diverges 
as £ ~ e J l T . 

At a normal first-order point, on the other hand, no 
critical behavior exists (i.e. all relevant L-independent 
length-scales remain finite), but the eigen-exponent d also 
shows up in finite-size corrections. In order to illustrate 
this, let us consider jllj a field-driven first-order tran- 
sition at h = (like e.g. in the two-dimensional Ising 
model below T c ). Assume that the central spin in a sys- 
tem of linear size L is pinned in the direction opposite to 
the external field h, and let us take for simplicity T ~ 0. 
If L is large enough, the bulk energy associated with 
h will be predominant, and m(L, h) will have the same 
sign as h. On the other hand if L is small, it will be 
energetically favorable not to break any bonds and thus 
all spins will point in the —h direction. The crossover 
length L* at which magnetization inversion happens can 
be estimated by equating the total field contribution L d h 
to the energy cost associated with the breaking a finite 
number of bonds around the pinned spin, and one obtains 
L* ~ h^ 1 ^ [jlT|] . So we see that there may be a diverging 
crossover size in a first-order transition, although there 
is no diverging correlation length in this case. 
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Thus the fact that finite-size corrections in the global 
observable I of small- world networks behave as L~ d is 
compatible both with a first-order transition and with a 
first-order critical point at p = 0. Clearly it is still neces- 
sary to prove or disprove the existence of critical behavior 
in small- world networks, and for this purpose we suggest 
to look at a bilocal observable (from which the analo- 
gous of a correlation function g(r) could be eventually 
defined) , instead of the global observable £. We study in 
this work the behavior of £(r), the average shortest-path 
distance between two points separated by an Euclidean 
distance r. 




shortcuts (dashed lines) have been connected between 



domly chosen points. The typical separation between nearby 
shortcut-ends is thus s ~ p~ 1//d . A disease starts to spread at 
t = from A with unit velocity, and for short times only the 
primary sphere (shaded) of radius t becomes infected. Each 
time a sphere hits a shortcut-end, a secondary sphere starts 
to grow from the other end of the shortcut, which is typically 
located at a distance L from there. 

The basic question we wish to answer is whether it is 
possible to identify a characteristic length £ (p) on small- 
world networks, such that: i) is a length scale dictat- 
ing the behavior of £(r), ii) £(p) only depends on p and 
not on system size L, when L — > oo, and, Hi) £ diverges 
at the critical point p c . 

Argollo et al [§) have recently given a simple asymp- 
totic argument to suggest that, when L — > oo, no finite 
characteristic length can be defined for bilocal observ- 
ables such as shortest-path lengths on small-world net- 
works. Newman and Watts Q on the other hand, re- 
cently claimed that the characteristic length for shortest- 
paths is the mean separation s(p) ~ p~ x l d between 



shortcut-ends. This quantity clearly satisfies criteria ii) 
and Hi) above, but we will show that it does not satisfy 
i), i.e. it is not relevant for shortest-path lengths on large 
systems. 

In order to test the ideas above, consider the spread 
of a disease in a (i-dimensional small-world system of lin- 
ear size L, as depicted in Fig. ED. For simplicity we work 
on the continuum. Assume that, at t — 0, an infectious 
disease starts to spread radially with constant velocity 
v = 1 from A. Let the sphere of radius t grown from 
A be called "primary sphere" of infected sites, or "indi- 
viduals". Because of the existence of shortcuts, which 
we assume are traversed in zero time, other (secondary) 
sources of infection will appear at times tk at random lo- 
cations in the available non-infected space. These times 
tk are the times at which an infection sphere (primary or 
not) hits one end of a shortcut. We call the spheres born 
at the other end of these shortcuts, "secondary spheres" . 
In this setting, the shortest-path distance £(A, x) from 
A to x is exactly the time t at which a point x is first 
infected. If a point at x is hit by the primary sphere 
first (before being hit by a secondary sphere), this will 
necessarily happen at time t = (1e{A.,x), where d(A,x) 
is the Euclidean distance from A to x. In this case the 
shortest-path distance £(A, x) is simply d E (A, x). If on 
the other hand a secondary sphere hits x at an earlier 
time then £(A, B) — t' < d E (A, x). 

In order to better uncover the physical meaning of a 
characteristic length r c for shortest-path distances £(r), 
let us consider a slightly modified disease-spreading prob- 
lem. Let us assume that the primary version of the dis- 
ease (the one that starts at A) is lethal, while secondary 
versions are modified when traversing shortcuts, and be- 
come of a milder, non-lethal form, which protects indi- 
viduals who are infected by it against the lethal form. 
Therefore noinfected individuals who are hit by the pri- 
mary sphere will be killed, while those being hit first by 
a secondary sphere become immune. Now assume that 
a characteristic length r c exists such that £(r) « r for 
r < r c and £(r) << r for r > r c . Because of the re- 
lationship between disease-spreading and shortest-path 
distances described above, one would conclude that a) a 
person at x is killed by the disease with a high proba- 
bility if d E {A, x) < r c , and b) a person at x is relatively 
safe if d E (A,x) > r c . Therefore r c , the characteristic 
length of shortest-path distances £(r), plays the role of 
a "safety distance" from the primary infection site A in 
this modified disease-spreading model. 

Newman and Watts recently provided a clever recur- 
sive equation for the total volume V(t) of infected sites 
as a function of time Q, and solved it exactly in one 
dimension. Their result reads 

Vit) = 1 s (>/* _ x) . (1) 
where s ~ is the average separation between ends of 
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different shortcuts. Notice that V(t) ~ t for t << s (only 
the primary sphere is important) while later V(t) ~ e*/ s 
when £ > s (proliferation of secondary spheres) . Because 
of this change in the behavior of V{t) at t w s, New- 
man and Watts M suggest that s is the equivalent of a 
correlation length. 

The total infected volume V(t) in ([!]) is the sum of the 
primary volume V\(t) plus the secondary volume ^(t). 
The primary volume is simply V\(t) = t d fld, where Q<i 
is the volume of a d-dimensional hypcrshpere of unit ra- 
dius. According to our discussion above, we are inter- 
ested in calculating the probability Ps(r) for an individ- 
ual located at a distance r from A to become infected 
by a secondary sphere before being hit by the primary 
sphere. Because the secondary volume V2 (t) is randomly 
distributed in the total available space L d , it is easy to 
calculate the probability P2(t) for an individudal to be- 
come infected by the secondary version of the disease 
at time t or earlier. This is simply the fraction of the 
total volume that is covered by the secondary volume: 
P2{t) = Vi{i)jL d 1 and does not depend on the Euclidean 
distance from the primary source at A. If the individ- 
ual does not become immunized by the secondary dis- 
ease first, it will certainly be killed by the primary in- 
fection at time t = d E {A,x), when the primary sphere 
reaches him. Thus the relevant quantity is P2{t) eval- 
uated at t = d E (A,x), which measures the probability 
Ps(d E (A,x)) for an individual at x to survive. In one 
dimension one obtains, using (0) and Vi{t) = 2t, 



P s (d E (A,x)) = — 1 < 



4d E (A,x)/s _ 



Ad E (A,x) 



- 1 



(2) 



which goes to zero when L — » 00 no matter how large 
d E (A,x) is. Thus there is no finite safety distance in 
the thermodynamic limit, meaning that r c diverges with 
system size. The size-dependence of the "safety distance" 
r c can be estimated from @ by equating P5 w constant, 
from which one obtains r c ~ sln(L/s). 

These heuristic arguments are confirmed by an exact 
calculation of £{r) Jl3|], as described in the following. 
Newman and Watts equation for V(t) || can be exactly 
solved in all dimensions, and using this solution it is pos- 
sible to calculate the average shortest-path distance £{r) 
between two points separated by an Euclidean distance 
r. The result is @: 

£ M f r for r < r c = p- 1 / d \og{p 1 / d L) ^ 
\ r c for r > r c ^ ' 

We confirmed this result numerically in one dimension, 
by measuring the average shortest path distance £(r) as 
a function of Euclidean separation r, for several values 
of L and p. Equation |Q) then means that the rescaled 
shortest-path distance £{r) = £{r) jr c should be a univer- 
sal function of the rescaled Euclidean distance f = r/r c . 



This is confirmed by the fact that all data displayed in 
Fig. H fall on a single curve. 

Notice that the characteristic length r c ~ s\og(L/s) is 
much larger than s, the mean separation between short- 
cuts. The reason why s is not a relevant length scale 
for shortest-path distances is not difficult to understand. 
Consider the shortest-path distance £(A, x) between A 
and some point x as in Fig. g. While the average dis- 
tance one has to travel from A in order to find a shortcut 
end is s, this shortcut has a typical length ~ L and thus 
will take us too far away from A to be useful in shortening 
the distance to x. 
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FIG. 2. Rescaled shortest path distance £(r)/r c ver- 
sus rescaled Euclidean distance r/r c , for one-dimensional 
small- world networks. The good collapse of these data onto a 
single curve clearly confirm Eq. (^) in the text. The charac- 
teristic length r c for shortest-paths behaves as p _1//d log(L d p) 
and thus diverges, for all p, with system size L. 

According to (^), only when d E {A, x) is larger than 
r c = s\og(L d p) is it more "convenient" to go from A 
to x through shortcuts than directly through the lattice. 
This result means that the typical shortest path in the 
regime of long distances r » r c contains 0(log(L d p)) 
shortcuts, since the typical "price" paid for each shortcut 
traversed is the average distance between shortcut ends, 
i.e. s. 

Thus s is not a relevant scale for shortest-path lengths 
£{r) and therefore it does not play the role of a correla- 
tion length as suggested recently J^]. A physically more 
appealing interpretation of s is to say B that it is a 
diverging timescale for the dynamic process of disease 
spreading. 

As might be clear by now, the impossibility to define 
a finite correlation length is a consequence of the fact 
that the locations of the secondary spheres are uncorre- 
cted with the location of the primary infection. In other 
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words, while the typical separation between ends of dif- 
ferent shortcuts s(p) is finite for L — > oo, the typical sep- 
aration between both ends of the same shortcut, which 
is an important quantity for shortest-paths, scales as L. 
As a consequence of this, the number of shortcuts with 
both ends inside a region of any finite size £ goes to zero 
when L —> oo, which already implies the impossibility 
to define a finite correlation length ||. Thus, although 
our demonstration here that no finite correlation length 
can be defined for shortest-path lengths is more rigor- 
ous, the same conclusion can be reached by using simple 
arguments H|. 

A different situation would certainly arise if shortcuts 
had a length-dependent distribution, for example if each 
site i is connected with probability p to another site j 
chosen with probability r~", where r is the distance be- 
tween i and j and a is a free parameter. For a — > 0, this 
model is the same as discussed here, while for a large 
enough one would only have short-range connections and 
thus there would be no logarithmic regime in £, even for 
p = 1. We speculate that the transition in p could be 
shifted to a finite p c for some intermediate a values. We 
are presently working on this extended model ]T^|. 
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